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Abstract 
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O \ A detailed study is made of four dimensional SU(2) gauge theory with static 

adjoint "quarks" in the context of string breaking. A tadpole-improved action 
is used to do simulations on lattices with coarse spatial spacings a^, allow- 
ing the static potential to be probed at large separations at a dramatically 
0^ , reduced computational cost. Highly anisotropic lattices are used, with fine 

temporal spacings at, in order to assess the behavior of the time-dependent 
effective potentials. The lattice spacings are determined from the potentials 
for quarks in the fundamental representation. Simulations of the Wilson loop 
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> 
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■ in the adjoint representation are done, and the energies of magnetic and elec- 



o 



trie "gluelumps" (adjoint quark-gluon bound states) are calculated, which set 



Q \ the energy scale for string breaking. Correlators of gauge-fixed static quark 

I propagators, without a connecting string of spatial links, are analyzed. Corre- 

lation functions of gluelump pairs are also considered; similar correlators have 
■ recently been proposed for observing string breaking in full QCD and other 

_ . models. A thorough discussion of the relevance of Wilson loops over other op- 

D , erators for studies of string breaking is presented, using the simulation results 



presented here to support a number of new arguments. 



^ ■ I. INTRODUCTION AND MOTIVATION 



Simulations of lattice QCD are increasingly dedicated to the goal of including the effects 
of sea quarks on many observables. One of the most distinctive signatures of sea quarks 
should be the elimination of the confining potential between widely separated valence quarks. 
Quenched simulations have demonstrated that color-electric field lines connecting a static 
quark and antiquark are squeezed into a narrow tube or "string" . In full QCD however the 
fiux-tube should be unstable against fission at large separations -R, where sea quarks should 
materialize from the vacuum and bind to the heavy quarks to form a pair of color-neutral 
bound states. 

It is perhaps surprising that some controversy persists in the literature as to whether 
"string breaking" has actually been observed in lattice simulations, and what techniques are 
required in order to convincingly demonstrate this phenomenon. This is despite extensive 
large scale simulations in unquenched QCD by several collaborations [0-§|. This problem 
has come under renewed attack in the last few years with several new viewpoints emerging 
as to the underlying cause of the difficulty in observing string breaking, and suggestions as 



1 



to the optimal approach for resolving this problem |17|-[T5|. These ideas have stimulated new 



work on string breaking in simulations of full QCD |[16Hl8i, and on a number of models that 
may shed light on string breaking . 



One suggestion to arise in the literature is that the Wilson loop operator, which has 
typically been used to study the static potential between heavy quarks, has a very small 
overlap with the true ground state of the system at large R, and hence is not suited to studies 
of string breaking P,[7|,pil[| . This has led to the consideration of other operators to study string 
breaking between static quarks, especially operators that explicitly generate light valence 
quarks at the positions of the heavy quarks |TT,T^. We note that another distinguishing 



feature of the operators proposed in Refs. [|n|)0 ^'^^^ they generate correlation functions 
which receive disconnected contributions, where a pair of heavy-light bound states propagate 
independently of one another. 



Another point of view was raised by one of us in Ref. where it was suggested that 
string breaking can indeed be seen using Wilson loops, but that it is essential to propagate 
the trial states over Euclidean times T of about 1 fm, the characteristic scale associated 
with hadronic binding. In contrast, typical studies of the static potential in unquenched 
QCD have been done on lattices with relatively "fine" spacings, limiting the propagation 
times to well under 1 fm, due to the high computational overhead. The use of coarse lattices 
enables much more efficient simulations of the static potential at the large scales relevant to 
string breaking. This was demonstrated in Ref. [|13| where an improved action was used to 
observe string breaking on coarse lattices, using Wilson loops in unquenched QCD in three 
dimensions, at dramatically reduced computational cost. 

In this paper we consider a number of the issues raised in the recent literature on string 
breaking, by doing simulations of quenched lattice QCD with valence "quarks" in the adjoint 
representation of the color group. There is a long history of lattice simulations of QCD with 
adjoint matter fields, which exhibits much of the physics of confinement of real QCD, and 



which also has a connection to supersymmetric physics [^. In particular an analogue of 
string breaking should occur in this model. For instance, the confining fiux-tube between 
a pair of heavy adjoint quarks should be unstable against fission at large R, where gluons 
can materialize from the sea to bind to the heavy quarks, forming a pair of color-neutral 
bound states dubbed "gluelumps" ||25|. Hence the potential Kdj(-R) between a pair of static 
adjoint quarks in quenched QCD should approach a constant at large R 

V;dj(i?^oo) = 2Mq„ (1) 

where Mgg is the energy of the lightest gluelump. Once again, despite much effort , 
string breaking using Wilson loops has not been seen in this model. However it has been 
suggested very recently that string breaking for adjoint quarks can be seen using correlators 
that explicitly generate valence gluons p0|-p2| . 



Here we undertake a detailed study of four dimensional SU(2) gauge theory with static 
adjoint quarks (this work was reported in unpublished form in Ref. |31[). We use a tadpole- 
improved gluon action to do simulations on lattices with coarse spatial spacings a^. We 
use highly anisotropic lattices, with fine "temporal" spacings at, in order to make a careful 
study of the time- dependent effective potentials. One lattice used here has = 0.36 fm and 
at = 0.10 fm, which provides an increase in computational efficiency of some two orders of 



magnitude compared to simulations of adjoint quarks that were done in Ref. [22] using an 
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unimproved action on lattices with spacings of about 0.1 fm. 

The lattice spacings are first determined from the potentials for quarks in the funda- 
mental representation. We then study the Wilson loop in the adjoint representation, and 
the masses of magnetic and electric gluelumps, which set the energy scale for string break- 
ing according to Eq. (|1|). We consider correlators of gauge- fixed static quark propagators, 
without the string of spatial links that is found in the Wilson loop; this is similar to correla- 
tors proposed in Refs. [^,0 as alternatives to the Wilson loop for observing string breaking. 
We also consider correlation functions for a pair of gluelumps, corresponding to states that 
explicitly contain valence gluons. 

We find that adjoint quark string breaking is extremely difficult to observe using Wilson 
loops, because of a very strong suppression of the signal, due to approximate Casimir scaling 
of the static potential. Despite the considerable computational advantage provided by using 
coarse lattices, we are limited to propagation times well below 1 fm. It is clear from our 
results that plateaus in the effective mass plots cannot be established at such short times. 
However it is also clear than a progressive "flattening" of the adjoint potential occurs as the 
propagation time is increased. 

By contrast saturation is readily observed in the effective potential naively constructed 
from the gluelump-pair correlator. Similar results for some of the observables studied here 



were recently obtained on "fine" lattices using unimproved actions [20-22]. The adjoint 



quark results in Refs. p0|^2| were interpreted as providing support for the picture p|,|7|JT0[| 
that Wilson loops are not suitable for studies of string breaking, while suggesting that string 
breaking can readily be observed using operators that explicitly generate states containing 
light valence particles . 

However we suggest instead that one must be more careful in defining the goal of "ob- 
serving string breaking" in Euclidean time correlation functions. In particular we note 
that operators which explicitly generate light valence quarks, such as proposed in Refs. 
]TT| , p]^P(]| -P^ , will automatically exhibit potentials that saturate, even in quenched QCD 
(here considering the theory with fundamental representation quarks). This is especially 
clear when the correlator has disconnected contributions. This can also be seen from re- 
sults presented in Ref. [|I9|, where correlation functions of exactly this form were studied 
in quenched QCD with fundamental representation quarks. If one goal of string breaking 
studies is to observe a distinctive feature of the effects of sea quarks, then clearly one must 
use an observable that distinguishes between the quenched and unquenched theories. 

We further argue that using an operator such as the Wilson loop (which might be thought 
of as creating "initial" states containing only heavy valence quarks) is well motivated by 
making an analogy with the process of hadronization, which is the physical process that 
is most often alluded to in the context of "string breaking." In hadronization the initial 
state consists of just two valence quarks at small separations, which then separate in real 
time, leading to the creation of light valence quarks at separations around 1-2 fm. The 
small overlap of the Wilson loop with the broken string state at much larger separations 
may be an irreducible feature of using Euclidean time correlators to mimic hadronization; 
a nonlocal state of two widely separated heavy quarks connected by a string is bound to 
have a small overlap with the state consisting of two widely separated heavy-light mesons. 
On the other hand, we find that the overlap of the Wilson loop with the broken string state 
is appreciable in a range of separations around the point at which string breaking actually 
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occurs. Simulations of Wilson loops at very large separations, where the overlap becomes 
very small, are not relevant to hadronization, since in this physical process the original 
quarks never get to such points with the string intact. 

If one considers the static potential for quark separations around 1-2 fm, which are 
physically motivated by the analogy with hadronization, then "string breaking" does indeed 
appear to be accessible in full QCD using Wilson loops. The important observation here is 
that propagation times of about 1 fm must be attained in order to resolve the broken string 



state |T3|. This can be more readily achieved by using improved actions to do simulations on 
coarse lattices; we estimate that an increase in computational efficiency of some two orders of 
magnitude over most recent studies in full QCD can reasonably be achieved. This viewpoint 
is also supported by the results presented here, given the behavior of the time dependent 
effective masses, even though actual string breaking in adjoint Wilson loops is not resolved. 
In this connection, we draw attention below to the somewhat surprising fact that string 
breaking in quenched QCD with adjoint quarks actually has a much higher computational 
burden than in unquenched QCD with real quarks. In our estimation all of the evidence 
supports the view that string breaking is accessible in large scale simulations of the Wilson 
loop in full QCD. 

The rest of this paper is organized as follows. In Sect. II we present the details of 
the improved gluon action, and the construction of the various correlation functions to be 
studied. The results of the simulations are presented in Sect. Ill, where each correlator is 
considered in turn. Finally, we present some further discussion and conclusions in Sect. IV. 



II. ACTION AND OBSERVABLES 



The tadpole-improved SU(2) action on anisotropic lattices used here was previously 
studied in Ref. following earlier work in SU(3) color 
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where P^j, is one half the trace of the 1x1 Wilson loop in the ix x v plane, R^v is one half 
the trace of the 2x1 rectangle in the x z/ plane, and where is the bare lattice anisotropy, 

V u-s / bare 

This action has rectangles Rgs' and Rst that extend at most one lattice spacing in the time 
direction. This has the advantage of eliminating a negative residue high energy pole in 
the gluon propagator that would be present if Rts rectangles were included. "Diagonal" 
correlation functions computed from this action thus decrease monotonically with time, 
which is very important for our purposes. The leading discretization errors in this action 
are thus of O(a^) and 0{asO?s). 

On an anisotropic lattice one has two mean fields Ut and Us- Here we define the mean 
fields using the measured values of the average plaquettes. Since the lattice spacings at in 



our simulations are small, we adopt the following prescription [p4 -|34| for the mean fields 
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Ut = 1, Us= {Pss'Y'^- (4) 

Observables in various representations of the gauge group can be easily computed from 
the measured values of the fundamental representation link variables using relations amongst 
the group characters. The Wilson loop Wj in the j-th representation is defined by 



where VjlUi] is the j-representation of the link t//, and L denotes the path of links in the 
Wilson loop. In the case of the adjoint Wilson loop of interest here, we have 

Wr{T,R) = \ (4 1 1^1/2 (T,i?) 1 2 - l) , (6) 

as can also be seen by using an explicit form for the adjoint representation matrices p5| 

Vl^\U] = ^Tiia'^Ua^U^), (7) 

and making use of the identity <jfj<Jki = '^{^u^jk — ^^ij^M) 



To enhance the signal-to-noise we make an analytical integration on time-like links 35 



/ d[U,]V,me-^' = '-fl^^AVi] I d[U,]e-^' (8) 

where kiVi denotes the sum of the 1x1 staples and 2x1 rectangles connected to the time-like 
link Ui [det{Vi) = 1]. This variance reduction was applied to time-like links for the Wilson 
loops and gluelump correlators. Equation (H) assumes that a given link appears linearly 
in the observable, hence we can only apply it to Wilson loops with R > 2, because of the 
rectangles Rst that appear in the improved action. 

An iterative fuzzing procedure was used to increase the overlap of the Wilson loop 
and gluelump operators with the lowest-lying states. Fuzzy link variables U^"'^ (x) at the nth 
step of the iteration were obtained from a linear combination of the link and surrounding 
staples from the previous step 

Ut\x) = ut'\x) + eY: f/f-'^(x)f/f"-^)(x + i)t//"-i)^a; + ?), (9) 

where i and j are purely spatial indices, and where the links were normalized to f/^f/ = / 
after each iteration. Operators were constructed by using the fuzzy spatial link variables in 
place of the original links. Typically the number of iterations n and the parameter e were 
chosen around (n, e) = (10,0.04) for Wilson loops and (n, e) = (4,0.1) for gluelumps. 

The gauge- invariant propagator G{T) of a gluelump can be constructed by coupling a 
static quark propagator Qt (product of temporal links) of time-extent T to spatial plaquettes 
Uq and Ut located at the temporal ends of the line 



G{T) = Tr {Uocx') Vf [Qt] Tr [u^a') . (10) 
Both magnetic and electric gluelump propagators were analyzed, by choosing appropriate 



linear combinations of spatial plaquettes at the ends of the static propagator |25|. For 
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the magnetic gluelump a sum of four plaquettes in a particular spatial plane is used, the 
sum being invariant under lattice rotations about an axis perpendicular the plane. For 
the electric gluelump a sum of eight plaquettes lying in two planes is used, the sum being 
invariant under rotations about an axis that lies in both planes. G{T) can be expressed in 
terms of fundamental representation link variables using Eq. (|^) 

G{T) = 2Tr (UoQtU^QI^) - Tr (Uo) Tr (^7^) . (11) 

We also study the correlation between two gluelump propagators, measuring the expec- 
tation value of the operator 

GGG(T,i?)=Gt(T;i?)G(T;0), (12) 

where static quark propagators Qt;o and Qt;R of time-extent T, and separated by a spatial 
distance R, are used in G(T; 0) and G(T; R), respectively. We also compute the off-diagonal 
entries in the gluelump pair- Wilson loop mixing matrix, given by the expectation value of 
the operators 

Ggw{T, R) = Tr (f/o;oa'^) Vf [Qt.o Tt.r QI,^] Tr {Uo.rct'') , (13) 

and 

Gwg{T, R) = Tr (Ut;rct') Vf Iq^ T^^ QtJ Tr (Ut;o(t') . (14) 



FojH and Tt-r are products of (fuzzy) links connecting the spatial sites of the heavy quarks at 
times zero and T respectively. The plaquettes Uq-q and Uq-r in the case of Ggw for example 
are connected to the ends of the static propagators at time zero. 

Finally, as an alternative to the Wilson loop, we compute correlators of gauge-fixed static 
quark propagators separated by a distance R, given by expectation values of the operator 

G'poty(T, R) = Tr {V,[Qt.,o]) Tr {V^[Qt;r\) . (15) 

This operator is similar to the Wilson loop in that it has only heavy quark propagators. 
It has been suggested ||,^ that this type of operator may have a larger overlap with the 
broken string state, since it does not have an explicit string of spatial links connecting 
the heavy quarks, unlike the Wilson loop. We measured Gpoiy in lattice Coulomb gauge, 
where J2^=i[Ui{x) — Ui{x — z)] = 0, which we implemented using an iterative steepest ascent 
algorithm with fast Fourier acceleration [R^ . 



III. RESULTS 

A. Lattice parameters and fundamental representation potentials 

Four lattices were studied in order to check the physical results for dependence on lattice 
spacing and input anisotropy. The four sets of simulation parameters are listed in Table |. 
Roughly 40,000 measurements were made for the observables on all lattices, skipping 10 
configurations between measurements (which results in very small autocorrelation times). 
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^0 


Us 


Volume 




as (fm) 


at (fm) 


AViVSas) 


0.848 


0.276 


0.7933 


10^ X 20 


1.02(1) 


0.361(8) 


0.102(2) 


0.078(2) 


0.848 


0.125 


0.8432 


8^ X 30 


1.17(1) 


0.494(20) 


0.072(3) 


0.154(2) 


0.600 


0.125 


0.7947 


8^ X 30 


1.14(2) 


0.606(40) 


0.086(6) 


0.220(2) 


0.500 


0.125 


0.7648 


8^ X 30 


1.12(1) 


0.689(30) 


0.096(4) 


0.231(2) 



TABLE I. Simulation parameters for the four lattices, and measured values of some lattice 
quantities. The bare anisotropies and the mean fields Us for tadpole improvement are shown 
(where ut = 1), along with the lattice volume in each case. Measured values of the lattice 
anisotropies ^ren are compared to the input anisotropies, as discussed in the text. Simulation 
results for the spatial and temporal spacings and at are given, as well as the relative errors AV 
in the off-axis potentials at R = y/Sas- 



We note that lattices with coarse spatial spacings were deliberately chosen, in an 
effort to probe the potentials at the longest physical quark separations possible, for the 
least computational cost. Lattice spacings around 0.4 fm have been studied by a number 
of authors (see e.g. Refs. []38| , |3^ , |32| ) . Using the lattice with as = 0.36 fm here provides 
an increase in computational efficiency of some two orders of magnitude compared to the 
simulations of adjoint quarks done on lattices with "fine" spacings in Ref. [Q. Still coarser 
spacings were also considered here and, although one would not necessarily advocate the 
use of such lattices for precision calculations, it is worthwhile to employ them in an effort 
to glean some information about string breaking, where much remains to be understood. 

We first present results for the fundamental representation Wilson loop, which are used 
to measure the renormalized lattice anisotropy ,^ren and to set the lattice spacing a^. These 
results demonstrate fairly good scaling and rotational symmetry restoration, as well as fairly 
small renormalizations of the input anisotropy, thanks to the tadpole improvement of the 
action. 

The renormalized anisotropy is determined by comparing the static potential atVxt, com- 
puted in units of at from Wilson loops W^t where the time axis is taken in the direction of 
small lattice spacings, with the potential a^V^y computed from Wilson loops W^y with both 
axes taken in the direction of large lattice spacings [p2|- ^ , ^ . The anisotropy is determined 
after an unphysical constant is removed from the potentials, by subtraction of the simulation 
results at two different radii 

^ _ atVxtiR2) - atVxtjRi) ^^g^ 

asVxy{R2) — asVxy{Rl) 

The anisotropies determined with Ri = \/2as and R2 = 2as are shown in Table |; results 
obtained with Ri = as are in excellent agreement with these estimates. The renormaliza- 
tion of the anisotropy is fairly small in all cases, especially as compared to the very large 
renormalizations for unimproved actions on lattices with comparable spacings . 

The spatial lattice spacing is then determined by fitting the fundamental representation 
potential to the form 

Vfit{R) = cTR-^ + c, (17) 
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FIG. 1. Fundamental representation potential for the action with = 0.49 fm: on-axis points 
(+), off-axis points (x). The dotted line shows the results of a fit of the on-axis points to Eq. (Ilj 
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FIG. 2. Time-dependent effective mass plot for the fundamental potential from the lattice with 
as = 0.49 fm and at = 0.072 fm. Each roughly horizontal line of points shows the effective mass at 
one separation R: on-axis points (-I-), off-axis points (x). 



8 




FIG. 3. Fundamental representation potentials in physical units from all four lattices: 
as = 0.36 fm (+), = 0.49 fm (x), = 0.61 fm (*), and = 0.69 fm (□). An additive 
renormalization in the energies has been adjusted so that the potentials agree at i? ~ 0.5 fm. The 
dotted line is the best fit to Eq. (^) for the lattice with = 0.36 fm. 
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FIG. 4. Effective mass plots for single electric (*) and magnetic (□) gluelumps. There are four 
pairs of plots corresponding to the four lattices, with the spatial spacing increasing from the 
bottom of the figure to the top. 
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FIG. 5. Gluelump physical mass splitting versus lattice spacing a^. The dashed curve shows a 
fit assuming O(a^) scaling violations. 



taking the physical value of the string tension to be -^/a = 0.44 GeV. The potentials were 
measured at on-axis separations, as well as at off-axis separations R/ag = V^, \/3, \/5, \/8, 
a/TS, vTS, and v^20- Symmetric combinations of the shortest spatial paths connecting two 
lattice points were used in the nonplanar Wilson loop calculations. Results for the lattice 
with ttg = 0.49 fm are given in Fig. |l], which shows fairly good rotational symmetry restora- 
tion, thanks to tadpole improvement. A quantitative measure of the symmetry breaking is 
obtained by comparing the simulation results for the potential with the interpolation to the 
on-axis data 



AV{R) = 



VsUR)-Vm{R) 
aR 



(18) 



Results for AV^ at i? = a/So,, for the four lattices are shown in Table |. Unimproved actions 
exhibit much larger rotational symmetry breaking effects on such coarse lattices ||32|| . 
Representative plots of the time-dependent effective potential 



y(T, i?) = - In 



\W{T-l,R) 



(19) 



are shown in Fig. |^. A reliable determination of the ground state potential in the fundamental 
representation can be made, with excellent plateaus in the effective mass plots, going out to 
propagation times near 1 fm, even at separations as large as 2.5 fm. 

The potentials from all four lattices are plotted together in physical units in Fig. |^. 
A Coulomb term is visible in this data, with fits to the on-axis data yielding coefficients h 
around 0.1 p^. Although the fundamental potentials on these coarse lattices are dominated 



by the confinement term, the results give some indication of fairly good scaling behavior. 
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B. Magnetic and electric gluelumps 



The effective mass plots for single electric and magnetic gluelumps exhibit good plateaus, 
as shown in Fig. ^. The electric gluelump is known to have the larger mass I^Sj. The 
gluelump energy Mqg is not physical as it must be additively renormalized due to the self 
energy of the heavy quark. However this renormalization should cancel in the difference 
between the electric and magnetic gluelump energies. A direct comparison of SMg^ with 
the static potential for a pair of adjoint quarks is also meaningful since the two quantities 
have equal self energies. 

Our results for the gluelump splittings on the four lattices are: 

{166 ±11 MeV, a, = 0.36fm, 

139 ± 15 MeV, a, = 0.49 fm, , . 

93 ± 16 MeV, a, = 0.61 fm, ^ > 

72 ± 18 MeV, a, = 0.69 fm, 

and are plotted versus lattice spacing in Fig. ^ For the sake of illustration, a fit assuming 
O(a^) scaling violations yields a continuum estimate of 

Melee " M^,g = (204 ± 16) McV, xVdof = 0.29, (21) 
while a fit assuming 0(a) scaling violations yields 

Melee - M^^^ = (273 ± 29) MeV, xV^of = 0.50. (22) 

These results are consistent with an estimate of the gluelump splitting in SU(2) color by 



Jorysz and Michael |25|, who found Meiec — Mmag = 203 ± 76 MeV, using a single lattice 



with a spacing of about 0.16 fm. 



C. Adjoint representation Wilson loops 

The determination of the ground state potential in the adjoint representation is much 
more difficult than in the fundamental case, due at least in part to a much larger overall 
energy scale in the adjoint channel. Effective mass plots for the adjoint potential on two 
lattices are shown in Figs. ^ and 0. Notice that the temporal spacing is smaller in Fig. |^, 
allowing one to more clearly see that plateaus in the effective masses at large separations 
have not been reached. The trend in the effective mass plots at large R is not inconsistent 
with the suggestion [|l3l that string breaking occurs at propagation times of about 1 fm. 



A typical procedure followed in the literature on string breaking is to estimate the ground 
state potential V^(-R) as equal to the effective potential V(T, R) at a small value of T, 
especially at large R, given the poor signal-to-noise in this region. The effect of choosing 
different fixed values of T for the determination of the potential can be seen by plotting 
V{T, R) versus R, for several choices of T. We show our data in this way for one lattice 
in Fig. H. There is a clear trend for the "potential" to flatten as T is increased, and this 
trend continues until the signal at large R is lost in the noise. The limitation to such small 
propagation times T ^ 0.4 fm at large separations introduces a significant systematic error 
in assessing whether the potential saturates. 
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FIG. 6. Adjoint Wilson loop effective mass plots for the lattice with Ug = 0.36 fm and 
Of = 0.10 fm. Plots are shown for several values of the on-axis quark separation, R = 1-3 (+), 
i? = 4 (□) and R = 5 (*), as well as at some off-axis points (x). The dashed lines show the la 
limits for twice the mass of the magnetic gluelump 2MQg on this lattice. 
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FIG. 7. Adjoint Wilson loop effective mass plots for the lattice with = 0.49 fm and 
at = 0.072 fm. Plots are shown for R = 1-3 {+), R = 4 (x), R = 5 {*). 
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FIG. 8. Adjoint Wilson loop static potential vs R at various fixed propagation times, for the 
lattice with = 0.49 fm and at = 0.07 fm: T/at = 2 (o), 3 (+), 4 (□), 5 (x), 6 (A). The dotted 
lines show la limits for 2MQg (magnetic). 



It is also useful to compare the fundamental potential with the adjoint one. In Fig. ^ 
we plot the two potentials in physical units from all four lattices, where we rescale the 
adjoint potential by 3/8, the ratio of SU(2) Casimirs for the two representations. There is 
good evidence for screening of the adjoint potential, compared to simple models of Casimir 



scaling ||25|, but again the extent of the screening (or possible saturation) cannot be reliably 



determined, due to the limitation to short propagation times. 



D. Gauge-fixed quark-antiquark correlator 



The correlation function between a pair of static adjoint quark propagators (cf. Eq. 
(p!5|)) was calculated in Coulomb gauge, in order to study a state without an explicit string 
of links connecting the heavy quarks. Similar correlators were suggested for observing string 
breaking in Refs. Results for the effective potential defined from Gpoiy(T, R) are shown 
for one lattice in Fig. [10|. The results obtained from this correlator agree well with the Wilson 
loop estimate of the potential, obtained at similar propagation times, giving neither a better 
nor a worse indication of string breaking. 



E. Gluelump-gluelump correlators 



Representative effective mass plots for the magnetic gluelump-gluelump correlator (cf. 
Eq. (|12D ) for one lattice are shown in Fig. |Tl]. At smaller separations the signal is clear 



but contains large excited state contributions; at larger separations the signal degrades, but 
the data at small T show more of a plateau. The resulting potentials from two lattices are 
compared in physical units in Fig. 
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FIG. 9. Comparison of the adjoint (*) and fundamental (A) potentials in physical units from 
all four lattices. The adjoint potential has been multiplied by a factor of 3/8 and shifted vertically 
to agree with the fundamental potential at i? 0.5 fm. The dashed lines show Icr limits for 2MQg 
(magnetic), after being rescaled and shifted vertically by the same amount as the adjoint potential 
(here using the results only at = 0.36 fm). 
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FIG. 10. Adjoint potential computed from gauge-fixed static quark propagators on the lattice 
with as = 0.36 fm and at = 0.10 fm. The dotted line shows 2MQg (magnetic), and the dashed 
curve shows the result of a fit to the potential computed from the adjoint Wilson loop on the same 
lattice and at comparable propagation times. 
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FIG. 11. Magnetic ghichimp-ghiclump effective masses for the lattice with ag = 0.36 fm and 
at = 0.10 fm, for two separations: R = 1 (lower points) and R = A (upper points). 




FIG. 12. Static potential estimates from the magnetic gluelump-pair correlator in physical units 
(with additive energy shifts): Qg = 0.36 fm (+) and Qg = 0.49 fm (x). 
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We also used a standard variational method p5[ to estimate the state of lowest energy 
in the 2x2 basis of states composed of a pair of heavy adjoint quarks connected to each 
other by a string of links (adjoint Wilson loop Wadj), and a pair of gluelumps. Consider the 
corresponding 2x2 correlation matrix Cij (cf. Eqs. (|12D-(|14D) 

(23) 



r (T m-( (^adj(T,i?)) {Ggw{t,R)) 



Representing the two basis states by |0j(i?)), the correlation matrix Cij(T,R) is written as 
a transfer matrix 

Q,{T,R) = {^,{R)\e-''^\<P,{R)). (24) 
We first find a linear combination |$(i?)) of basis states 

mR)) = J2a,{R)\<P,{R)) (25) 

i 

which maximizes 

iHRmR)) ■ ^^^^ 

This requires the solution of the eigenvalue problem 

Cij{T\ R)aj{R) - X{T\ R)Cij{0, R)aj{R) = 0. (27) 



We choose to optimize the variational state by solving Eq. (|27D at a small time T*, otherwise 
numerical instabilities may arise due to large statistical errors in Cij{T* , R), especially at 
large R. 

We then evolve the correlation function Eq. (^) for the variational state to a larger time 
T, in order to filter out our final estimate of the ground state energy A(T, R) = Q-^oi^)'^ _ 
The overlaps cf{R) = {(f)i{R)\0{R))^ / {(j)i{R)\(j)i{R)) of the basis states |0i(i?)) on the ground 
state \0{R)) can be estimated according to 

at sufficiently large T (Eq. (^) at finite T provides an upper bound on the overlaps). 

The results of this diagonalization procedure are as might be expected. Figure |T3| shows 
the estimate of the ground state potential in physical units from two lattices. The estimated 
overlaps of the Wilson loop and gluelump-pair states with the variational estimate of the 



ground state are shown in Fig. |14[ There is a rapid cross-over in the ground state as 
determined in this basis, from the Wilson loop at smaller R to the gluelump-pair state at 
larger R. The rapid cross-over suggests that the gluelump-pair state becomes dominated by 
disconnected contributions when the energies of the states near 2MQg. 
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FIG. 13. Variational estimate of the ground state energy in physical units for two lattices (after 
additive shifts in the energies): Ug = 0.36 fm (+) and = 0.49 fm (x). The trial state was typically 
determined at T* = at, which was then propagated to a time T ~ 4at to obtain the results shown 
in the figure. Also shown are the la lines for 2MQg (magnetic). 
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FIG. 14. Diagonalization results for the lattice with Ug = 0.49 fm, showing the overlaps of the 
two basis states with the ground state as functions of R, estimated from Eq. (p^): Wilson loop 
(+) and gluelump-pair (*). 
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IV. SUMMARY AND FURTHER DISCUSSION 



In this paper we made a through analysis of adjoint quark physics in the context of string 
breaking. Three trial states were investigated as candidates for observing string breaking: 
the adjoint Wilson loop, a pair of gluelump propagators, and a pair of gauge-fixed static 
quark propagators. The fundamental representation potentials were used to measure the 
lattice parameters, and electric and magnetic gluelump masses were calculated in order to 
set the energy scale for string breaking. A number of techniques were used to maximize 
the efficiency of the calculations, including fuzzing, variance reduction, and fast Fourier 
accelerated gauge fixing. Particularly important was the use of coarse, highly anisotropic 
lattices from tadpole-improved actions. The lattice with = 0.36 fm and at = 0.10 fm for 
example gives an improvement in computational efficiency of some two orders of magnitude 
compared to simulations of adjoint quarks done in Ref. on lattices with spacings of about 
0.1 fm. Large lattice anisotropics provided more data points for analysis of the Euclidean 
time evolution of correlation functions. This proved to be especially important in analyzing 
adjoint Wilson loops for even moderate physical values of R, due to the rapid decay of the 
signal. 

Diagonalization in a two state basis, corresponding to adjoint Wilson loops and the 
gluelump-pair operators, reveals a static potential which saturates at 2MQg near 1.5 fm. 
Similar string breaking distances have been suggested for full QCD and have been 



observed in other theories, including three-dimensional QCD with dynamical fermions []T3 
At small quark separations the potential rises linearly, with a slope of about | of the fun- 
damental potential, consistent with Casimir scaling. Saturation in the potential obtained 
from the diagonalization occurs over a very small range of separations R. These results 
are in qualitative agreement with recent calculations done on fine lattices using unimproved 
actions p0|-[2^. In addition, results obtained here using correlations between gauge-fixed 
static quark propagators agree well with the Wilson loop estimate of the potential, giving 
neither a better nor a worse indication of string breaking. 

As discussed in Sect. I these results, taken at face value, might be interpreted as evidence 
that Wilson loops are not suitable for studies of string breaking, and that string breaking 
can be readily observed using operators that explicitly contain fields for the light particles. 
However we raised several cautionary observations about this viewpoint. In particular, we 
pointed out in Sect. I that operators which explicitly generate a pair of light valence parti- 
cles automatically exhibit potentials that saturate. In the case of QCD with fundamental 
representation quarks, this means that operators of this type would demonstrate "string 
breaking" even in the quenched theory. This is particularly clear when the correlation func- 
tion receives disconnected contributions. In the case of the gluelump-pair correlator (cf. Eq. 
(P^)), we have 

{0\G\T; R)G{T- 0)|0) = (0|G'(T; 0)|0)2 + . . . , (29) 

where |0) is the vacuum state, and where the ellipsis denotes the contributions of non- 
vacuum insertions between the two operators in the correlator. Hence one is guaranteed to 
find an effective mass of 2MQg at modest R using this correlation function. This also helps 
to explain the rapid cross-over that is observed in variational calculations using Wilson loops 
and such states |lTT| , p^ , pO| -p^ , when the energies of the states approach the broken string 
energy. 
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Exactly the same behavior must occur if operators that exphcitly generate hght valence 
quarks are used in simulations of QCD with fundamental representation quarks. This can 
also be seen from some recent work using such operators in quenched QCD Thus 
operators of the type that have recently been proposed in Refs. [0,0 for observing string 
breaking will yield qualitatively the same behavior in both quenched and unquenched QCD. 

This raises the important question of exactly how one defines the goal of "observing 
string breaking." Perhaps one can advocate two points of view. From one point of view, 
one would say that correlation functions give spectral information only: the ground state 
energy versus quark separation. In this view such information is useful, for example, in 
characterizing the effects of quenching on meson- meson interaction energies ||4T| , |i2[] , but not 
in relation to processes such as hadronization, which can only occur in the presence of 
dynamical fermions. Then one might say that the difficulty in seeing "string breaking" with 
the Wilson loop at very large separations indicates nothing more than that this operator 
has a poor overlap with the ground state in this regime. If spectral information is all that 
one is interested in, then one might not be too concerned that a particular operator shows 
qualitatively similar behaviour in the quenched and unquenched theories. 

A second point of view, which we raised in Sect. I, is that one can use certain correlation 
functions to make an analogy with hadronization. One might argue that this may be done 
by considering operators that do not generate light valence particles in the trial state. In 
particular, the operator should exhibit "string breaking" only in unquenched QCD (here 
considering the theory with fundamental representation quarks). The Wilson loop is such an 
operator, while operators which explicitly generate light valence particles are not. Moreover, 
if a goal of string breaking studies is to observe a distinctive feature of the effects of sea 
quarks, then one should only consider observables that distinguish between the quenched 
and unquenched theories. In a sense the first point of view discussed above is concerned with 
"static" spectral information, that is qualitatively similar in the quenched and unquenched 
theories, while the second point of view makes contact with a "dynamical" process that is 
unique to the unquenched theory. 

In making contact with hadronization one is interested in using the Wilson loop to 
measure the potential for separations R not much larger than the distance at which the 
string breaks, i? ^ 1.5 fm, since in the actual physical process the original quarks never get 
to larger separations with the string intact. In this region the overlap of the Wilson loop 
with the ground state appears to be appreciable, judging from Fig. |14|, and from results 
presented in Ref. [0]. However one must still push the calculation to propagation times 
of about 1 fm, in order to properly resolve the broken string. This is where the difficulty 
actually lies in observing string breaking using Wilson loops. 

To date most simulations of full QCD have been done on lattices with relatively fine 
spacings, making it computationally very challenging to reach the length scales i? ^ 1.5 fm 
and propagation times T ^ 1 fm relevant to string breaking. The use of coarse lattices with 
improved actions allows a much more efficient probe of this regime, as was recently demon- 
strated by one of us in Ref. [jT^, where this approach was used to resolve string breaking 
with Wilson loops in unquenched QCD in three dimensions. An increase in computational 
efficiency of some two orders of magnitude is possible using lattices with spacings between 
0.3 fm and 0.4 fm, compared to most simulations that have been done so far in unquenched 
QCD. Some work in this direction has recently been reported in Ref. ]T7|. 
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Unfortunately calculations of the adjoint Wilson loop in this string breaking regime did 
not prove to be feasible in this paper, with propagation times limited to well below 1 fm, 
even with coarse lattices. On the other hand it is clear that adjoint Wilson loops exhibit a 
potential that progressively "flattens" at longer propagation times. Moreover, the trend in 
the effective mass plots from the adjoint Wilson loop at large R is not inconsistent with the 
suggestion that string breaking would be observed at propagation times of about 1 fm. 

In this context it is interesting to estimate the size of the adjoint Wilson loop signal 
relative to the fundamental one, and to compare the computational cost of these simulations 
to those of unquenched QCD. If we assume roughly Casimir scaling of the potential just below 
the string breaking distance, then the ratio VFadj/W^fund of the adjoint to the fundamental 
Wilson loops in SU(3) color goes like 

l^adj/mund ~ exp [-(I - l)aRT] ^ 10-^ (30) 

using ^/a = 0.44 GeV for the fundamental representation string tension, and R ~ 1.5 fm 
and T ^ 1 fm for the scales relevant to string breaking (the ratio is yet smaller, by an 
order of magnitude, in SU(2) color). This is to be compared with the roughly two orders of 
magnitude increase in the cost of simulating dynamical quarks compared to quenched sim- 
ulations. (Note that the gluelump-gluelump correlator p0|^2| does not show a comparable 
suppression of the signal, which is entirely a consequence of not removing the disconnected 
contributions from the correlator). Hence, while the adjoint representation is interesting 
as a probe of confinement and supersymmetric physics, it is not a cost effective means of 
mimicking hadronization in full QCD. 

On the other hand the results presented here do lend support to the general picture 
that Wilson loops should in fact exhibit "string breaking," as an analogy to hadronization. 
Moreover string breaking should be accessible in real QCD with the computational power 
currently available in large scale simulation environments, especially if coarse lattices with 
improved actions are used. 
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